CONFORMAL INVARIANCE IN RANDOM CLUSTER MODELS. 
I. HOLOMORPHIC FERMIONS IN THE ISING MODEL. 



STANISLAV SMIRNOV 



1. Introduction 



It is widely believed that many planar lattice models at the critical temperature are conformally 
invariant in the scaling limit. In particular, the Ising model is often cited as a classical example of 
conformal invariance which is used in deriving many of its properties. 

To the best of our knowledge no mathematical proof of this assertion has ever been given. 
Moreover, most of the physics arguments concern rectangular domains only (like a plane or a strip) 
or do not take boundary conditions into account. Thus they give (often unrigorous) justification 
only of the Mobius invariance of the scaling limit, arguably a much weaker property than full 
conformal invariance. Of course, success of conformal field theory methods in describing the Ising 
model provides some evidence for the conformal invariance, but it does not offer an explanation or 
a proof of the latter. 

It seems that ours is the first paper, where actual conformal invariance of some observables for 
the Ising model at criticality (in domains with appropriate boundary conditions) is established. Our 
methods are different from those employed before, and allow us to obtain sharper versions of some 
of the known results. Moreover they allow the construction of conformally invariant observables 
in domains with complicated boundary conditions and on Riemann surfaces. Ultimately we will 
construct conformally invariant scaling limits of interfaces (random cluster boundaries) and identify 
them with Schramm's SLE curves and related loop ensembles. These extensions will be discussed 
in the sequels jl5l [TB] . Though one can argue whether the scaling limits of interfaces in the Ising 
model are of physical relevance, their identification opens possibility for computation of correlation 
functions and other objects of interest in physics. 

We consider the Fortuin-Kasteleyn random cluster representation of the Ising model on the 
square lattice ST? at the critical temperature. This representation, briefly reviewed below, studies 
random clusters, which are clusters of the critical percolation performed on the Ising spin clusters 
at the critical temperature. The spin correlations can be easily related to connectivity properties in 
the new model. Every conflguration can be described by a collection of interfaces (between random 
clusters and dual random clusters) which are disjoint loops that fill all the edges of the medial 
lattice. 

As a conformally invariant observable we construct a "discrete holomorphic fermion". In a 
simply connected domain Q with two boundary points a and b we introduce Dobrushin boundary 
conditions, which enforce the existence (besides many loops) of an interface running from a to 6, 
see Figure [TJ We show that the expectation that this interface passes through a point z taken with 
fermionic weight (i.e. a passage in the same direction but with a 27r twist has a relative weight — 1, 
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whereas a passage in the opposite direction with a counterclockwise vr twist has a relative weight 
—i - see Figure [3]) is a discrete holomorphic function of z. Moreover, as the step of the lattice 
goes to zero, this expectation, when appropriately normalized, converge to a conformally covariant 
scaling limit, namely V^, where $ is the conformal map of to a horizontal strip. 

The approach is set up for a random cluster model with a general value of the parameter q £ [0, 4] , 
and a parafermion observable of spin a = 1 — ;| arccos(y^/2), conjecturally converging to {^'Y 
in the general case. The Ising case corresponds to g = 2 and a = 1/2. Besides a priori estimates 
(which are well-known in the Ising case), we make essential use of the Ising-specific properties in 
two places: to establish discrete analyticity of an observable, and to show that being a solution 
of the discrete Riemann boundary value problem, it converges to its continuum counterpart. For 
the latter we see possibilities for a proof in the general case, albeit more difficult. So it seems 
that the only essential obstacle to proving conformal invariance of all random cluster models lies 
in establishing discrete analyticity of the observable concerned. For the Ising case this is done by 
proving discrete analogues of Cauchy-Riemann relations, where partial results can be obtained for 
all random cluster models. 

The two sequels [151 EB] are concerned with the construction (on the basis of one observable) 
of conformally invariant scaling limits of one interface and full collection of interfaces respectively. 
They are, more or less, applicable to all random cluster models for which conclusions of this first 
part, in particular Theorem 12.21 can be established. In the Ising case the law of one interface 
converges to that of the Schramm-Lowner Evolution with n = 16/3. 

These results were announced and the proofs were sketched in [T3], where one can find some of 
the ideas leading to our approach. Another notable case when this approach (or rather a parallel 
one) works is the usual spin representation of the Ising model at critical temperature, leading to a 
similar observable (related to a conformal map to a halfplane), and to Schramm-Lowner Evolution 
with K, = 3. On a rectangular lattice, exactly the same notion of discrete analyticity arises. This 
will be discussed in a separate paper. 

Similar observables were constructed before by Kenyon [5] and by the author \13\ I12j . The 
work of Kenyon concerned dimers on the square lattice (domino tilings) and so by the Temperley 
bijection gave conformally invariant observables for the Uniform Spanning Tree (corresponding to 
the random cluster model with q = 0) and Loop Erased Random Walk. Kenyon's considerations 
are close in spirit to ours, in fact repeating his constructions for the Fisher lattice one is led to 
similar observables. Since SLE was not available at the moment, the identification of interfaces 
had to wait till the work ^ of Lawler, Schramm and Werner. Nevertheless Kenyon was able to 
rigorously determinate several exponents and dimensions, and some of his results go beyond the 
reach of SLE machinery. We constructed |13t [T2] conformally invariant observables for the critical 
site percolation on the triangular lattice which also bear some similarity to ones in the current 
paper. Unfortunately that proof is restricted to the triangular lattice, so the question of conformal 
invariance remains open for the percolation on the square lattice (which corresponds to the random 
cluster model with q = 1). 

The paper is organized as follows. In Section [2] we state our theorem. We start the proof by 
introducing a new notion of discrete analyticity in Section [3l and then show that it is satisfied by 
an Ising model observable, which we construct in Section HI Finally in Section [5] we show that the 
discrete observable has a conformally covariant scaling limit. In the proof we use an a priori estimate 
for the Ising model, which follows from (a weak form) of known magnetization estimates; this is 
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Figure 1. Loop representation of the random cluster Ising model. Weight of the 
configuration is proportional to (y^)* ^"°p^^ with q = 2. The sites of the original 
Ising lattice are colored in black, while the sites of the dual lattice are colored in 
white. Loops separate clusters from dual clusters, which are also pictured, the 
former in bold. Under Dobrushin boundary conditions besides a number of loops 
there is an interface running from a to 6, which is drawn in bold. 



discussed in Appendix |Al Some of the more technical results about discrete harmonic functions 
are reviewed in Appendix [Bl 
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2. Statement of results 

We work with the Fortuin-Kasteleyn random cluster model with a particular emphasis on the 
critical Ising case, corresponding to parameter values q = 2 and p = \/2/(\/2 + !)• For a general 
introduction to the Ising and random cluster models consult the books [H El [9]. 

The random cluster measure on a graph (a simply connected domain on the square lattice in 
our case) is a probability measure on edge configurations (when each edge is declared either open 
or closed), such that the probability of a configuration is proportional to 

p# open edges ^-j^ p)^ closed edges g# clusters 

where clusters are maximal subgraphs connected by open edges. The two parameters are edge- 
weight p £ [0, 1] and cluster- weight q G [0, oo), with q £ [0, 4] being of interest to us. For a square 
lattice (or in general any planar graph) to every configuration one can prescribe a random cluster 
configuration on the dual graph, such that every open edge is intersected by a dual closed edge and 
vice versa. See Figure [1] for two dual configurations with open edges pictured. It turns out that 
the probability of a dual configuration becomes proportional to 

# dual open edges /-i \#: dual closed edges # dual clusters 

P* P*) q 1 

with the dual to p value = p*{p) satisfying p*/(l — p*) = ^(1 —p)/p- For p = p^d := y^/ (■v/9 + 1) 
the dual value coincides with the original one: one gets psd = ipsdj* and so the model is self-dual. It 
is conjectured that this is also the critical value of p, which was only proved for g = 1 (percolation), 
q = 2 (Ising) and q > 25.72. For these and other properties of the random cluster models consult 
Grimmett's monograph [3j. 

We will work with the loop representation, which is perhaps the easiest to visualize. The cluster 
configurations can be represented as (Hamiltonian) loop configurations on the medial lattice (a 
square lattice which has as vertices edge centers of the original lattice), with loops representing 
interfaces between cluster and dual clusters - see Figure [H It is well-known that for p = p^d the 
probability of a configuration is proportional to 

(^)# loops ^ 

with q = 2 m. the Ising case. 

We introduce Dobrushin boundary conditions: wired on the counterclockwise arc ha (meaning 
that all edges along the arc are open) and dual-wired on the counterclockwise arc ah (meaning that 
all dual edges along the arc are open, or equivalently all primal edges orthogonal to the arc are 
closed) - see Figure [H For the loop representation this reduces to introducing two vertices with 
odd number of edges: a source a and a sink h. Then besides a number of closed loop interfaces 
there is a unique interface running from a to 6, which separates the cluster containing the arc ah 
from the dual cluster containing the arc ba. See Figures [U H] for typical configurations. 

Note that Dobrushin boundary conditions are usually formulated for the spin Ising model and 
amount to setting plus and minus spin boundary conditions on two arcs correspondingly, thus 
creating an interfaces between two spin clusters. Since we need an interface between two random 
clusters, we formulated similar conditions in the random cluster setting. Our version of Dobrushin 
boundary conditions is equivalent to setting plus boundary conditions on one arc and free on the 
other in the spin setting. 

The model makes sense for q = Q as well, and is equivalent to the uniform spanning tree model. 
Indeed, setting g = prohibits loops, so we consider configurations containing only an interface 



Conformal invariance in random cluster models. I. 



5 



from a to b, which are weighted uniformly since the number of open edges is always the same. Those 
configurations are easily seen to be equivalent to spanning trees on the original spin lattice, rooted 
on the arc ab. We are mostly interested in the Ising case of q = 2 and p = psd = \/2/(\/2 + 1). 
Though it is known that this value of p is critical, we won't use it in the proof. For other values of 
q our proof works to large extent, also for the self dual value p = Psd, and in principle one can try 
to use this in establishing its criticality. 

Define spin by cr := 1 — ;| arccos(y^/2). Note that for the Ising case a = ^. Let F[z) be the 
expectation that the interface 7 passes through a point z taken with a complex weight: 

F{z) := E(x^g^(a;) • exp{-iaw{j,b ^ z))) . 

Here w denotes the winding or (the total turn) of 7 from b to z, measured in radians. For the Ising 
case an additional 2tt turn of the curve before reaching z changes the weight by a factor of —1, 
see Figure [3l The formula above gives F at the edge centers (of the medial lattice, where the loop 
representation is defined), and we extend it to all of O in a standard piecewise constant way. Exact 
definition can be found below. 

We start by proving in Section [5] the following 

Proposition 2.1. For the Ising model in a given lattice domain the function F{z) is discrete 
holomorphic and satisfies a discrete analogue of the Riemann Boundary Value Problem (with a = ■^J 

(1) Im {F{z) tangent {z f) = 0. 

The continuum problem is solved by / = {^')'^ , where ^ is the conformal map of i7 to a horizontal 
strip, with a and b mapped to the ends. In our case normalization will produce a trip of width 2. 
After some technicalities in Section [5] we show that F converges to its continuum counterpart: 

Theorem 2.2. Suppose that as the lattice mesh 6j goes to zero, the discrete domains 0,j on the 
lattices Sjl? (with points aj, bj on the boundary) converge to the domain $7 (with points a, b on 
the boundary) in the Caratheodory sense. Then for the Ising model the corresponding functions 
Fj = F{z,Q.j,aj,bj,5jl?) converge uniformly away from the boundary: 

(2) 8rF,^f = {^'y , 

where a = 1/2. 

Remark 2.3. Caratheodory convergence is defined as convergence of normalized Riemann uni- 
formization maps on compact subsets. 

Namely we fix a point w £ Q, and let (j) (or (j)j) be conformal maps from the unit disk D to ^7 (or 
Q-j) such that points 0, 1, C (or 0, 1, Q) are mapped to w, a,b (or w, aj, bj) or corresponding prime 
ends. We say that 0,j converge to 0, if converge to (p inside D and Cj tends to (. 

It is easy to see that Hausdorff convergence of the boundaries implies Caratheodory convergence 
and that solution to the Riemann boundary value problem (CP for z inside ^l, being defined in the 
terms of Riemann maps, is uniformly continuous as a function of Q in the topology of Caratheodory 
convergence. 

A variation of our proof seems to work for g = as well, and most of it can be worked out for 
other values of q, though sometimes in a different way. Hopefully the missing part of the discrete 
analyticity statement will be worked out someday: 

Conjecture 2.4. Proposition \2.1\ (with an appropriate, possibly approximate, discrete analyticity) 
and Theorem \2.2\ hold for all values of q £ [0,4]. 
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3. Discrete analyticity revisited 



We will identify lines through the origin with unit vectors (complex numbers) belonging to them. 
For a line I or equivalently a vector a G £ we denote by Proj (-F, €) = Proj (-F, a) the orthogonal 
projection of a complex number F on the line i. Note that for a unit vector a 



Consider the square lattice bl? (possibly rotated). By a lattice domain £l we mean some collection 
of vertices joint by edges such that all vertices have even number of edges. In our application we 
will also allow half-edges (usually two, their middle ends will become the source a and the sink 6). 

By distance between two points (when speaking of moduli of continuity of functions, etc.) inside 

we will mean the distance in the inner metric. 

If for some vertex all four edges are present, we call it an interior vertex^ while vertices with two 
edges we call boundary vertices. If for a square all four vertices are interior, we call it an interior 
square. To the lattice domain Q we associate a planar domain Q which is the union of all interior 
squares. We will assume that those domains are connected and simply connected. 

Color the lattice squares in a chessboard fashion. We orient every edge e, turning it into a unit 
vector (or a complex number) e with the orientation chosen so that the white square is on the left 
and the black one on the right. Then to the edge e prescribe a line l{e) in the complex plane which 
passes through the origin and the square root of the complex conjugate of the vector e, considered 
as a complex number (note that the choice of the square root is not important). 

Without loss of generality we can assume that the lattice edges are parallel to the coordinate 
axis (otherwise all lines are rotated by a fixed angle). Then horizontal edges correspond to the 
lines with argument (defined up to tt) or 7r/2 (in the chessboard order), whereas vertical edges 
correspond to tt/A or 37r/4. See Figure [2j 

A given vertex v has 4 neighboring edges. If we go around v counterclockwise, the line cor- 
responding to the neighbor with each step is rotated counterclockwise by 7r/4 (so the full turn 
corresponds to a rotation by vr, which preserves the line but reverses directions). 

Definition 3.1. We say that a function F defined on vertices is preholomorphic or discrete analytic 
in a domain ft if for every edge e £ Q orthogonal projections of the values of F at its endpoints on 
the line £{e) coincide. We will denote this common projection by F{e). 

Remark 3.2. In the complex plane holomorphic (i.e. having a complex derivative) and analytic 
(i.e. admitting a power series expansion) functions are the same, so the terms are often inter- 
changed. Though the term discrete analytic is in wide use, in discrete setting there are no power 
expansions, so it would be more appropriate to speak of discrete holomorphic ( or discrete regular ) 
functions. We prefer the term preholomorphic, which was common at one point, but seems to have 
gone out of use. 

Remark 3.3. The more commonly used discrete analyticity condition asks for the discrete version 
of the Cauchy-Riemann equations 



to be satisfied. Namely, for a lattice square with corner vertices NW , NE, SE, SW (starting from 
the upper left and going clockwise) one has 



(3) 




diaF = idaF 



(4) 



F{NW) - F{SE) = i {F{NE) - F{SW)) 
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It is easy to check that our property implies the more common one on vertices, but does not fol- 
low from it. Moreover, our property is equivalent to the more common property for the function 
restricted to horizontal edges. 

Definition 3.4. We will say that a preholomorphic function F solves the Riemann Boundary Value 
Problem ([I]) in the domain if for every boundary vertex v with two edges, projections of F{v) on 
the lines corresponding to these edges have the same modulus. 

Remark 3.5. Indeed, since projections of F{v) on lines corresponding to two edges coincide, F{v) 
belongs to the bisector of those lines. Equivalently the value of F at every boundary vertex v is 
parallel to the reciprocal of the square root of the tangent vector t(v) (or rather a discrete approx- 
imation - the vector orthogonal to the bisector of the angle between two edges from v). This is a 
discrete analogue of the Riemann Boundary Value Problem (OP and our main goal will be to show 
that solutions to the discrete problem in the limit solve the continuum one. 

We will solve the problem ^ by "integrating" the square of F, which is not so easy since is 
no longer preholomorphic. 

Lemma 3.6. Let F be preholomorphic in a domain 0,. Then up to a constant there is a unique 
function H = Hp defined on the lattice squares in or adjacent to $7 and such that for any two 
adjacent squares, say black B and white W separated by the edge e, one has 

(5) H{B)-H{W) = |F(e)|2 . 

In applications we will chose the constant so that H is zero on a square immediately below 6. 
Note that H is defined on the dual lattice to 5J?, which is the usual case for discrete derivatives 
or primitives. 

Remark 3.7. Values of the argument of F on edges are such that the function 25H is a discrete 
analogue of the indefinite integral Im J F'^dz, which we will establish in Appendix{^ equation II 33]): 
if u and v are centers of two squares with a common corner z, then 

2 6 [H{v) - H{u)) = Im {F(zf {v - u)) . 

Proof. It is sufficient to check that when one goes around an interior vertex, increments of H add 
up to zero. Suppose that the edge neighbors of the vertex v are E, S, W, N in clockwise order 
starting from the right. Then the sum of increments when we go around v is 

(6) ± \FiE)\^ T \Fm' ± \HW)\' T , 

with signs depending on the choice of chessboard coloring. By construction i{E) _L i{W), i{N) _L 
i{S). Since F{E), F{W), F{N), F{S) are orthogonal projections of F{z) on the corresponding 
lines, by Pythagoras theorem 

\Fiv)\^ = |F(F)|2 + \F{W)\^ = |F(iV)|2 + . 

Thus the sum of the increments is equal to 

±mE)\' + \F{w)\' - \Fm' - \Hs)\') = ±mv)\' - i^wp) = o , 



and indeed vanishes. 



□ 
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Figure 2. An example of a lattice domain O with boundary conditions creating 
an interface from a to b, which is drawn in bold. The lattice squares are colored in 
the chessboard fashion, with black corresponding to the sites of the original Ising 
lattice and white to the sites of its dual. Near the edges we write the arguments of 
the corresponding lines. Note that running from a to 6 the interface always follows 
the arrows and has black squares on the left. 



Denote H restricted to black squares by H'\ We define the discrete Laplacian by 

AH'{B) := J2 {H{Bj) - H{B)) , 

j 

where the sum is taken over four black squares Bj - neighbors of the black square B, touching it at 
vertices. Similarly we define the Laplacian for the restriction H'^ to white squares. We say that a 
function is discrete (sub/super) harmonic if its Laplacian vanishes (is non- negative/non-positive). 

Lemma 3.8. // F is preholomorphic, then for an interior white square W with corner vertices 
NW, NE, SE, SW (starting from the upper left and going clockwise) we have 

(7) AH'"{W) = -\F{NE)- F{SW)\'^ = - \F{NW) - F{SE)\'^ <0, 
so is superharmonic. Similarly, for an interior black square B 

(8) AH\B) = \F{NE) - F{SW)f = \F{NW) - F{SE)f >0, 
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SO H is subharmonic. 

Remark 3.9. It is clear that subharmonicity on black squares is equivalent to superharmonicity on 
white ones. Indeed, the definition of H is such that it is always increased when we pass from white 
squares to black. If we reverse the colors, we will arrive at the function —H and subharmonicity 
will become superharmonicity. 

Proof. A computation leading to equation (I7|8p is possible, since increment of H across some vertex, 
say NW, can be written in terms of the projections of F{NW) on various lines, and so ultimately 
in terms of F{NW). So we can express the Laplacian AH{B) in terms of the values of F at 
four neighboring vertices. But projections of F{NW) on two lines corresponding to upper and 
left edges of the square coincide with those of F{NE) and F{SW) correspondingly. So we can 
express F(NW) and similarly F{SE) in terms of F{NE) and F{SW). The resulting formula for 
the Laplacian is quite simple. 

The computation is rather lengthy, so we present it in the Appendix O However there are 
several arguments why we should arrive at a simple result. Since we deal with squares of absolute 
values of projections, we arrive at some real quadratic form in F{NE), F{NE), F{SW), F{SW). 
Symmetries of our setup imply that this form is invariant under the rotation by vr which yields the 
exchange F(NW) <-> —F{SE), and under change of F by an additive constant (this follows e.g. 
from the equations ()29l30l31l30p ). Such a form is necessarily given by const • \F{NE) — F{SW)\ . 
Note that any value of the constant would do, leading to a sub- or super- or harmonic function. □ 

By boundary arcs we mean parts of which are not separated by "distinguished" points (i.e. 
ends of half-edges). In our usual setup there are two boundary arcs, ah and ba (with points given 
counterclockwise). By values of Hp on the boundary we mean its values on the outside squares 
adjacent to fi. 

Lemma 3.10. // a preholomorphic function F solves the problem then Hp is constant on the 
boundary arcs. 

Proof. Go along a boundary arc over the squares adjacent to the domain. Let B and B' be the 
centers of two consecutive ones (say of black color), they touch at a vertex v, and are separated 
from 17 by the edges e and e' emanating from v. Then 

H{B) - H{B') = |F(e)|2 - \F{e')\^ = , 

So H is indeed constant along the arc. □ 

4. Discrete holomorphic spin structure 

We consider loop representation of the random cluster Ising model at critical temperature. The 
discrete domain $7 on the lattice 51? is as discussed above, with a "source" a and a "sink" b. Thus 
for every configuration w besides loops there is a curve 7 = ^{ijo) joining a and 6, which we will call 
the interface, see Figure [H Rotate the lattice in such a way so that an only edge incoming into h 
from Q. points to the right. 

We round the corners of the loops so that there are no sharp turns, see Figure [21 The loops can 
be connected at a vertex in two different ways, like near the vertex v in Figure O which are more 
clearly distinguished if we draw rounded loops. Note that the interface can pass through a vertex 
twice, utilizing two rounded corners - see the left part of Figure O Color the lattice squares in 
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weight 1 —i —1 

Figure 3. Values of the complex weight W {'y,z) for different passages of the in- 
terface 7 through z. 

a chessboard way, so that standing at b and facing the domain Q, we have a black square on the 
right and a white one on the left. The black squares correspond to the sites of the original Ising 
lattice, while the white squares correspond to the sites of the dual one. An interface going from a 
to b always has black squares on the left and white on the right, so it can arrive to a point z only 
from one direction (and not from the opposite one). Thus we can prescribe to every point a vector 
which is tangent to all interfaces passing through it from a to b. Since it has black square on the 
right, for points on the edge e it coincides with the vector e discussed above. 

Recall that to every point z (centers of edges and rounded corners are important) we prescribe 
a line i{z) in the complex plane which passes through the origin and the square root of the corre- 
sponding vector (the choice of the square root is irrelevant). For edges this agrees with the scheme 
discussed above, see Figure [21 

A given vertex v has 8 neighboring corner or edge centers. If we go around v counterclockwise, 
the line corresponding to the neighbor with each step is rotated clockwise by vr/S (so the full turn 
corresponds to a rotation by vr, which preserves the line but reverses directions). 

For two points z, z' on an interface 7 we will denote by w[^z z') = w{'y, z — > z') the winding 
(i.e. the total turn) of the curve 7 as it goes from z to z' , measured in radians. 

For an interface 7 we define the complex weight W at point z G 7 by 

W{'y, z) := exp (^(7> a ^ z) + w{'j, b ^ z) — w{j, a — > b)) 
Note that ^(7, a ^ z) — w{'y, b ^ z) = ^(7, a ^ b) and so 
(9) Ty(7,2;) = exp ( — u;(7, 6 ^ z) 



2 

The values of the complex weight are illustrated in Figure [3] 

Lemma 4.1. For a point z and any realization of the interface 7 the complex weight W{'^,z) 
belongs to the line i{z). 

Proof. When the interface is traced starting from b the property is easily checked by induction. At 
the center of the first edge the complex weight is equal to 1, and so belongs to the line through 1. 
When the interface turns by ±9, the winding w{'y, 6 — > z) is increased by ±9, whereas ^(7, a ^ z) 
is decreased by ^9. So the complex weight changes by a factor of exp(— z/4- {— ^ 9 ziz 9)) = 
exp {^i9/2). On the other hand, the line £{e) is also rotated by ^9/2 since it passes through the 
complex conjugate of the square root of the corresponding tangent vector which is rotated by ±9. 
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Here we use that the interface when traced from b always has black squares on the right and so 
goes in the opposite direction to the tangent vector. □ 

We will work with points z which are either "centers of corner turns" (near every vertex there 
are 4 such points - see Figured]) or centers of edges. For corner and edge points we can write 

Here n{z z') = 71,(7, ^ z') is the number of it^ turns with sign the curve 7 makes going from 
z to z' and A = exp(— ivr/S). Note that for corner points w{a z) differs from n{a — > z) ■ tt/2 by 
ib7r/4 (the last half-turn before reaching z) but the difference enters w{a ^ z) and w{b ^ z) with 
opposite signs and so cancels out. 

Remark 4.2. As was mentioned before, the choice of weight is such that a relative weight of the 
interface with an additional 2tt twist around z is —1. Indeed, such a twist forces each of the two 
halves a ^ z and b ^ z to make four tt /2-turns, so the weight for one Tr/2-turn is X = exp(— zvrS), 
which satisfies = —1. 

Remark 4.3. Taking A = exp(z7r/8) instead, one arrives at discrete anti-analytic functions. 

Define a function F at all "centers of corner turns" (near every vertex there are 4 such points - 
see picture) by 

F(c) := E (xeg^H • W{j{u;),c)) • 2 cos ^ . 

Similarly define F for all centers of edges by 

F(e) := E(xe6-,H-W^(7(w),e)) . 

Different normalization arises because there are more corners than edges per vertex. 

Remark 4.4. The given definition of F for edge centers works well only for the square lattice at 
criticality (which is perhaps the most interesting case). As an alternative one can start with our 
definition for corner centers, and use the equation [1^) to define F for edge centers. This approach 
gives the same function in our setting, but also generalizes to non-critical values of p and to other 
lattices. 

With corners rounded, the interface can go through a vertex f in 4 different ways, passing through 
one of the 4 neighboring corners Cj . For an interior vertex v we define F as 

F{v) := 5;F(c,)/2. 

j 

One can rephrase this as saying that 

(11) F{v) = E (x„g^(^) • W(-/{io),v)) • cos ^ , 

where all passages of the interface through v (there might be up to two) are counted separately. 

Lemma 4.5. For an interior vertex v the values of F at its 8 neighbors are orthogonal projections 
of F(v) on 8 corresponding lines. 

Remark 4.6. A boundary vertex v has only two neighboring edges, say e and e' . We define F(v) 
as a unique complex number which has orthogonal projections F{e) and F{e') on the corresponding 
lines i{e) and i{e'). It follows that F is preholomorphic in 
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Remark 4.7. The proof uses that we have a square lattice at v, hut the only global information 
needed is that the graph is planar. So if we define F on some planar graph which has square lattice 
pieces, it will he preholomorphic there. In a sequel we will discuss generalizations of preholomorphic 
functions to general planar graphs. 

Proof. When going around v clockwise the hne is rotated by vr/S with each step, thus hnes cor- 
responding to antipodal neighbors are at angle of 47r/8 = 7r/2, and so are orthogonal. Hence the 
values of F at two antipodes are orthogonal, and are orthogonal projections of their sum on the 
corresponding lines. So we will be proving in fact a stronger property, namely that 

(12) F{NW) + F{SE) = F{NE) + F{SW) = F{W) + F{E) = F{N) + F{S) = F{v) , 

where in each of the pairs two terms are orthogonal. Here starting from the right and going 
clockwise we denote 8 neighbors of v by E, SE, 5, SW, W, NW, N, NE. 
Recall that by definition 

F{v) = {F{NE)+F{NW) + F{SW) + F{SE)) /2 , 

so to establish the identity (I12p and the Lemma it is sufficient to show that the sum of values of F 
at two antipodal neighbors is the same for 4 such pairs of antipodes. 

Define an involution a; i— > on loop configurations, which results from the rearrangement of 
connections at the point v. For the random cluster formulation it corresponds to opening/closing 
the edge going through v. To prove the linear identity it is sufficient to show that each pair uj,uj' 
of configurations makes identical contributions to all 4 "antipodal" sums. 

Consider some pair of configurations, say oj and uj' . If the curve 7(w) does not pass through v, 
neither does 7(1-1-''); and all contributions are zeroes. 

Otherwise both curves pass through v. Trace either of the curves from a until it first arrives to 
the neighborhood of v. Since it has black squares on the left it can arrive from 2 possible directions, 
similarly when traced from b can arrive from 2 other directions. 

There are 4 possible topological pictures for the arrivals, but they are all analogous, so we will 
work out one of them. Assume that the half starting from a arrives from the west, while the half 
starting from b arrives from the south (such picture is possible for a half of the vertices v, for others 
the curve traced from a would arrive from a vertical direction). In one of the curves, say 7(0;), the 
two traced halves are immediately joined (and there is also a cycle passing near v), whereas in the 
other, 7(w'), this cycle is included into the curve. See Figures U] and [5l Then out of corner points 
7(0-1) contributes only to F{SW), say a term X (weight of all cycles times the complex weight). 
On the other hand the curve 'y{uj') out of the corner points contributes to F(^NW) and F{SE) 
only. The contributions differ from X by a factor of l/\/2, since the number of cycles decreased 
by one. Moreover, the phase changes, since compared to 7(1-1-') reaching SW the curve 7(0;') winds 
by additional +tt/2 when reaching NW (coming from the half originating in b) and by additional 
— 7r/2 when reaching SE (coming from the half originating in a). Correspondingly the complex 
weights change by factors of and A^, see equation (fTO]) . So 7(0;') contributes to F{NW) + F{SE) 
a term 

X-{\^ + X^)/V2 = X . (e-^i2 + e^f 2) = X • 2 cos (^) /V2 = X . 

So 7(0;) contributes X to the second sum in (fT2]) . while not contributing to the first, whereas 'y{uj') 
contributes X to the first sum in (fT2]) . while not contributing to the second. We conclude that the 
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Figure 4. Rearrangement at a vertex v: we only change connections inside a small 
circle marking v. Either interface does not pass through v in both configurations, 
or it passes in ways similar to the pictured above. On the left the interface (in bold) 
passes through v twice, on the right (after the rearrangement) it passes once, but a 
new loop through v appears. The loops not passing through v remain the same, so 
the weights of configurations differ by a factor of ^/q = \/2 because of the additional 
loop on the right. To get some linear relation on values of F, it is enough to check 
that any pair of such configurations makes equal contributions to two sides of the 
relation. 



configuration lo' 



configuration to 



to b 



to a Nw j^NE ^ 



to a W_ NW^\NE ^ 



SW\*SE 



to b 



Figure 5. Schematic drawing, representing the change in the interface after the 
rearrangement at a vertex v. Going from a and 6 to u the interface might make 
some number of turns, which won't influence our reasoning, since it changes the 
weight of both configurations by the same factor. Note that since a and b are on 
the boundary, for topological reasons the interface can go from N to E only on one 
side of V. 
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first two sums coincide: 

F{NW) + F{SE) = F{NE) + F{SW) . 
The same (but messier) reasoning shows that the two remaining sums share the same value. Per- 
haps this is best summarized by the following table which shows contributions of two configurations 
to the values of F at various neighbors of v: 





NW 


SE 


NE 


sw 


N 


S 


W 


E 













X 





X\ 


x\ 





2cos(7r/8) 


2cos(7r/8) 


J 




xx^ 










X\ 


XX 


xx-^ 


V2 


V2 


2cos(7r/8)V2 


2cos(7r/8)V2 


2cos(7r/8)V2 


2cos(7r/8)\/2 



Using that A = exp (— 'i7r/8), an exercise in trigonometry one checks that numbers in 2 x 2 
squares bordered by double lines always sum up to X. Thus taken together a; and to' make identical 
contributions to all 4 antipodal sums in (jl2p . 

Alternatively we can finish the proof by deducing that values of F on the neighboring edges 
are also projections of F{v). To that effect we write the value of F at the northern edge in 
terms of northwest and northeast corners. Consider some edge e emanating from v with the 
corresponding line ^{e) passing through a vector a. Let c and c' be the two adjacent corner 
points. The corresponding lines are rotations of £(e) by ±7r/8, passing through vectors Aa and Aa 
corr esp ondingly. 

Note that the interface passes through e if and only if it passes through exactly one of the points 
c and c'. Taking into account the difference in complex weight and normalization, and using the 
formula ([3]) for projections, we write 

(13) F{e) = {\F{c) + AF(c')) / (2 cos(^/8)) 

= (AProj (F(t;),AQ) + AProj {F{v),\a)) / (2 cos(7r/8)) 

= (A {F{v) + {XafFiv)) + A [F{v) + {\afF{v))) / (4cos(7r/8)) 

= {\F{v) + WF{v) + \F{v) + \a^F{v)) / (4cos(^/8)) 

= {F{v) + a^F{v)) (A + A) / (2(A + A)) = Proj {F{v), a) , 

thus finishing the proof. □ 

Lemma 4.8. For every positive r there is a function 5r{x) : M+ M+ such that lim^-^^o ^rix) = 
and if two neighboring squares B and W are r away from at least one of the boundary arcs ah or 
ha, then 

\H{B)-H{W)\ < 6r{6) . 

Remark 4.9. Here 6 is the lattice step. Note that the only way the shape of enters into the 
estimate is via r. The Lemma essentially means that the restrictions of the function H to black 
and white squares are uniformly close to each other when we are away from a and b. 

Remark 4.10. The Lemma is derived from the fact that F ^ away from a and 6 as 5 — > 0. 

Since preholomorphic F is uniquely determined by its boundary conditions, there should be an Lsing- 
independent proof, using only discrete analyticity and boundary conditions. Unfortunately, we were 
not able to find a simple one. 
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Proof. If an edge e separates the squares B and W, by definition 
(14) \H{B) - H{W)\ = |F(e)p < P (e E 7)^ 

so we can take as our function 6 the square of the similar function in Lemma lA.li □ 



Lemma 4.11. The function F satisfies the Riemann Boundary Value Problem ([I]). Moreover 
H = on the (counterclockwise) boundary arc ah and H = 1 on the (counterclockwise) boundary 
arc ha. 

Proof. Let v he a, boundary vertex with incoming edges e and e'. Then all possible interfaces 7 
arrive at e from a (or from b) with the same winding, so W{^{uj), e) is independent of uj. Therefore 

|F(e)| = |E (xe6-,H • Wi-f{io),e)) \ = \W{-f{Lo),e) E {xeejH) | = P (e G 7(^)) • 

Similarly |-F(e')| = P (e' G 7('-i^))- But since there are only two edges out of v, an interface passes 
through e if and only if it passes through e', so 

|F(e)| =P(eG 7(w)) =P(e' €7(0;)) = |F(e')| , 

and F satisfies the Riemann Boundary Value Problem ([T|), proving the first statement of the Lemma. 

By Lemma 13.101 it follows that the function H is constant on the boundary arcs. Let u and v be 
the centers of squares immediately below and above b. Recall that we chose H (which is determined 
up to an additive constant) so that H{u) = 0. Thus ff = on the (counterclockwise) boundary 
arc ab. Every interface passes through b, and furthermore has the same complex weight at b. So 

H{v) = H{u) + |F(6)|2 = + |E (X6g^ • W{^, b))\' = |1^(7, b) E (xte^)!' = P (6 G 7(^))' = 1 • 

Therefore H = \ on the (counterclockwise) boundary arc ba. □ 

When establishing the SLE connection in the sequel p!5], we will need the "martingale property" 
of F with respect to the interface. Consider the interface 7 as a random curve drawn from a 
to b with some parameterization. Let t < s be two stopping times (actually we do not need a 
full parameterization at this point). For the time t (and similarly for s) we denote by 7(i) the 
corresponding curve point, and by 7[0,t] the part of the curve from 7(0) := a to 7(t). When 
speaking of domain SI \ 7[0,t], we will actually mean its component of connectivity containing b. 

Then the following holds 

Lemma 4.12. Let z be a lattice vertex such that for every realization of the interface z G $7\7[0, s]. 
Then for every realization 0/7(0, t] 

F{z, n \ 7[0, t], 7(t), b) = E^[t,,] {F{z, Q \ 7[0, s], 7(5), 6)|7[0, t]). 

The conditional expectation above is taken over all possible continuations j[t, s] of the interface 
until the time s assuming the part 7[0,t] is given. 

Proof. With a fixed 7[0,t] consider the remaining part 7' of the curve 7 in the domain Q \ 7[0,4]. 
Plugging into the definition (jlip of F in this domain the formula ([9]) for the complex weight we 
obtain 

F{z,n\j[0,t],j{t),b) = E(X) , 
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where the random variable X is given by 

--w{j',b^ z) j • cos - . 
Now the desired formula is the law of total expectation applied to X and the curve j[t, s]: 



□ 



5. Passing to a limit 



In this Section we prove the Theorem 12. 2i To derive convergence we will use only discrete 
analyticity and boundary values of F, and appeal directly to the Ising connection only in quoting 
Lemma 14.81 As discussed in Remark 14.101 the latter should have an Ising-independent proof. So 
essentially the Theorem 12.21 can be restated as a Theorem about preholomorphic functions solving 
the Riemann boundary value problem ([1]). 

We work with a sequence of lattice domains, which approximate a given domain 0,. Consider 
a sequence of lattice domains with distinguished points aj, bj and with lattice steps 5j. Let 
Fj = F{z, Qj,aj, bj, bjl?^ be the expectation as defined above and denote := Hp.. 

Assume that 6j —>■ and Qj,aj,bj ^^—^ 0,a,6 as j — > oo. We use Caratheodory convergence 
of domains, which is the convergence of normalized Riemann uniformization maps on compact 
subsets. 

Recall that $ is a mapping of O to a horizontal strip M x [0, 2] , such that a and b are mapped 
to =Fc>o. Note that since <I> is uniquely defined up to a real additive constant, its derivative, and 
hence the right hand side in ([2]) are uniquely determined. Recall Remark 12.31 that solution to the 
continuum Riemann boundary value problem ([T]), the function is Caratheodory-stable. 

Remark 5.1. The convergence (0j in the Theorem \2. S\ holds on the boundary of wherever it is 
a horizontal or vertical segment. Since complex weight on the boundary is independent of configu- 
ration, we conclude that for such z € dVl 



from which one deduces that random cluster intersects the (smooth) boundary on a set of dimension 
1 /2 and that for the spin Ising model at criticality magnetization on the boundary is proportional 
to (with a specific factor). 

We start by establishing convergence of H^s: 

Lemma 5.2. Away from a and b there is a uniform convergence: 

(15) Hp^ ^ h:= Im$ . 

Proof. Remove the union V of some neighborhoods of a and b, then there is a positive r such that 
remaining parts of the boundary arcs ab and ba are at least 4r-apart. Then all points in $7 \ y 
are at least 2r away from at least one of the arcs ab and ba. We conclude that because of the 
Caratheodory convergence, for sufficiently large j all points in in Qj \ V are at least r away from 
at least one of the (discrete) arcs ab and ba. 



Conformal invariance in random cluster models. I. 



17 



Hence by Lemma 14.81 we have uniform convergence 

(16) sup \ H {B j) - H{Wj)\ =: 5j ^-^0, 

for neighboring squares Bj, Wj G Vlj \ V . 

Considering Hj restricted to black and white squares we obtain functions i/j" and HJ (super- 
harmonic and subharmonic correspondingly by Lemma 13. 8p . If we extend from black {HJ 
from white) lattice vertices to whole of $7 in any reasonable way (e.g. making constant on lattice 
squares), then ()16p can be rewritten as convergence in the uniform norm on Q \ 



(17) 



Tjb ttw 



=: 6, '^0 

n\v,oo 



Let Hj be a discrete harmonic function on black squares with boundary values given by h, define 
Hj' similarly. Then 

(18) + 25 j > + 5j > HJ - 6j > Hj - 25 j 

on the boundary, and hence inside domain i7 since the four functions involved are superharmonic, 
harmonic, harmonic and subharmonic correspondingly. Together with (|17|) this means that 



n\v,oo 

and since Hj^s converge to h by the Lemma lB.3| so do -ff/s. □ 

If Hj^s are harmonic, the Theorem immediately follows. Indeed, derivatives of Hj admit an 
integral representation (in terms of Hj itself), so uniform convergence of Hj implies uniform con- 
vergence inside Q of ^Hj and hence its square root, i.e. Fj. For a general approximately harmonic 
Hj this doesn't work, but in our case we can use that appropriate restriction of Hj is subharmonic 
and that Fj is exactly preholomorphic. 

First we will need the following compactness estimate: 



Lemma 5.3. Let U be a subdomain compactly contained in ^l, and denote by Uj its discretizations 

2 



with mesh 6j. Then 



is uniformly bounded. 

Remark 5.4. Note that the expression above is essentially norm of Fj/y^. 

Proof. Note that when we jump diagonally over a vertex z, the function Hj changes by ReFj or 
ImFj depending on the direction. It follows that it is enough to prove uniform boundedness of 

d,^\VH,\ 

where V denotes discrete difference gradient of Hj restricted to black or white vertices. 

From now on, we will work with Hj on the "black" sublattice, having centers of black squares 
as vertices. In particular, A will denote the corresponding Laplacian. Recall that restriction of Hj 
to this lattice is subharmonic, i.e. AHj > 0. 
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Fix a square Q of side length I such that a nine times bigger square is contained in Vt. 
We win denote by the same letter the lattice approximation Q = Qj of Q, whose side will have 
L = Lj := l/6j edges. It is enough to find for every such Q a uniform in j bound on 

(19) SjY,\^H,\. 

Qj 

Recall that on 9Q one has uniform convergence Hj^H. Denote by Hj the discrete harmonic 
function on 9Qj having the same boundary values as Hj. Then H*^H on d9Q, and by the 
Lemma ED iVi7*=^ViJ on Q. Summing over Q we infer that 

Qj 

converges to the area integral of \'VH\ and in particular is bounded. Thus to bound (|19|) it is 
enough to bound 

Since Hj and H* have the same limit, their difference converges uniformly to zero: 

(20) sup \H* — Hj \ = Oj —>■ when j oo. 

9Q 

Denote by Gj{x, y) = Ggg^. (x, y) the discrete Green's function for the square 9Qj, with lS.xG{x, y) = 
AyG{x,y) = 5x-y and vanishing on the boundary of 9Q x 9Q. Note that it is negative inside 9Q. 

Using discrete analogue of the Riesz representation formula we can write a subharmonic function 
Hj — Hj as a convolution of its Laplacian (which coincides with that of Hj) with the Green's 
function: 

(21) Hj{x) - H*{x) = AHj{y)G{x,y). 

J/69Q 

Taking difference gradient, we arrive at 

(22) V {H,{x) - H*{x)) = AH,{y)VxG{x,y). 

S/G9Q 

Using the Lemma lB.41 we write: 



Y\^{H,{x)-H;ix))\ = Y: 



Y AH,{y)VxG{x,y) 

yG9Q 



ye9Q x&Q 



^ const ^ A/7,(,) Y |G(x,y)| = ^ E E ^H,{y)G{x^y) 
ye9Q xeQ xeQyedQ 

const v-^ \ TT / \ rr* / N I ^ const / Oj 
— — — \-Hj\x) — Hj [x) I < const L Oj = 



J 



proving the Lemma. □ 
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Lemma 15.31 implies (by Theorem V.12.a in ^ applied to the primitives of -Fj's ) that the se- 
quence is precompact in the uniform topology on any compact subset of 17. To show its 
convergence to / it is sufficient to establish convergence to / of a uniformly converging (say to ^f) 
subsequence. Uniform convergence implies convergence of the (discrete) integral of Fj/^j to that 
of g'^. The imaginary part of the former is given by Hj + const, which also converges to Im$, so 
we conclude that the two limits are equal: 

Im$ = lim Hj = Im J + const. 

Since the function involved are analytic, equality of imaginary parts implies that they are equal up 
to a constant. Differentiating and taking the square root, we conclude that g = ^/W, thus proving 
the Theorem 12.21 

Appendix A. A priori estimates 

We use an estimate on the modulus of continuity of our function F, which essentially states 
that the interface cannot be space filling. It reduces to rather weak (compared to what is known) 
magnetization estimates, which ultimately can be retrieved from the (old) literature. However 
it seems that "modern" proofs are elusive and would be worthwhile. It is also possible to build 
everything on the basis of discrete analyticity, without appealing to properties of the Ising model. 
But for now we present a sketch of the proof using classical results and assuming unlike in the rest 
of the paper knowledge of the basic properties and techniques of the Ising model. 

Lemma A.l. For every positive r there is a function 5r{x) : M+ — > M-|- such that lim^-^o Sr{x) = 
and if an edge e is r away from at least one of the boundary arcs ah or ha, then 

P(e e 7) < 6r{5) . 

Proof. Denote by B and W the centers of two neighboring squares (black and white), separated 
by an edge e. If an edge e belongs to the interface, then B is connected by a cluster to the arc ha 
and ~ by a dual cluster to the arc ah (since the interface separates the two). Assume that the 
edge e is at least r away from the boundary arcs 6o, the case when it is away from ah is treated 
similarly with clusters replaced by dual clusters (which leads to the same result since the model is 
self-dual). Denote by Q the square with side length r/2 centered around e, by our assumption it 
does not intersect the boundary arc ha. Then by monotonicity we can write 

P(e G 7) = 

=P {B connected by a cluster to the arc ha and W connected by a dual cluster to the arc ah) 
<P {B connected by a cluster to the arc ha) 

<P {B connected by a cluster to dQ inside 17 with ha wired and ah dual — wired) 
<P {B connected by a cluster to dQ inside with haU ah\Q wired and abCiQ dual — wired) 
=P {B connected by a cluster to dQ inside 0, with dO, \ Q wired and d^l D Q dual — wired) 
<P {B connected by a cluster to dQ inside 0, with 17 n dQ wired and dO, n Q dual — wired) 
<P {B connected by a cluster to dQ inside Q with dQ wired ) 

= magnetization at B in the Ising spin model in Q with " + " boundary conditions on dQ. 
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The right hand side, the magnetization at criticahty, was computed by Kaufman-Onsager and 
Yang [HITl], and for a fixed r it tends to zero with the mesh 6 tending to zero, proving the Lemma. 

Note that magnetization was computed to tend to zero hke 5^^^, but we do not need this stronger 
statement. On the other hand, after convergence to SLE is estabhshed, one can obtain even better 
asymptotics for the magnetization, deriving a conformally covariant constant in front of 6^^^. □ 

Appendix B. Estimates of discrete harmonic functions 

In this Section we collect the needed facts about discrete harmonic functions. Let (5jZ^ be a 
sequence of lattices with mesh 6j tending to zero. For a domain U we denote by Uj the corresponding 
lattice domain in SjT^'^. As usual, functions on a lattice domain are thought of as defined in the 
whole domain, say by piecewise constant or linear continuation. 

A classical fact says that solution of the discrete Dirichlet problem converges uniformly away 
from the boundary to the solution of the continuum one. The following Lemma can be found in 
the seminal paper [2] by Courant, Friedrichs and Lewy. It can also be deduced from the random 
walk interpretation of harmonic functions (which is also discussed in [2]). 

Lemma B.l. Let 0,j be a sequence of lattice approximations to a smooth domain with the mesh 
6j tending to zero. Let {hj} be a sequence of discrete harmonic functions on lattice domains flj and 
h be a harmonic function on Q, with continuous boundary values. If hj converge uniformly to h on 
dQ, then inside U away from the boundary hj and its partial discrete derivatives (i.e. normalized 
differences) are equicontinuous and converge uniformly to their continuum counterparts, i.e. h and 
its partial derivatives. 

In our case the lattice domains approximate in the Caratheodory (rather than in the Hausdorff ) 
sense, furthermore the boundary values are discontinuous. We will deduce a suitable result using 
the following well-known estimate: 

Lemma B.2. There exist an increasing positive function e on M4. with limx.^o+ e(x) = 0, such 
that the following holds. Let H be a non-negative bounded discrete harmonic function in a simply- 
connected domain O with boundary values equal to zero on Q D B{z,r) and at most one elsewhere. 
Ifdist{z,dn) < S, then H{z) < e{5/r). 

This is a weaker version of discrete Beurling's estimate e{5/r) = const y^^/r. It can be refor- 
mulated in terms of the hitting probabilities for the simple random walk and is found in Kesten's 

Now we can prove the needed version of the convergence Lemma: 

Lemma B.3. Suppose that as the lattice mesh 5j goes to zero, the discrete domains Vtj (with 
points aj, bj on the boundary) converge to the domain $7 (with points a, b on the boundary) in the 
Caratheodory sense. Let hj be a discrete harmonic function on 0,j with boundary values on the 
arc ajbj and 1 on the arc bjaj. Then inside discrete functions hj converge uniformly to their 
continuum counterpart h, which is harmonic in $7 with boundary values on the arc ab and 1 on 
the arc ba. 

Proof. Being harmonic with bounded boundary values functions hj for large j are equicontinuous 
inside 17 by [2j or [17j - see inequality ([23|) below. Thus it is enough to show that any subsequential 
limit, say h', coincides with h. 
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Fix small R > such that two balls B{a, R) and B(b, R) are disjoint, and denote their union by 
W. Let r < 2R be the distance between the arcs ab\W and ba \ W. Take 5 < r/2 and let il.^ be 
a subdomain of O with smooth boundary which is (5/2-close to the boundary of fi. Let and b^ 
be two points on which are 5/2 close to a and b correspondingly. Caratheodory convergence 
implies that for large enough j subdomain is contained in and its boundary dQ^ is contained 
in the (5- neighborhood of the boundary of (The opposite inclusion might fail if Qj contains long 
fjords of fixed diameter, which however disappear in the Caratheodory limit if their width tends to 
zero). 

Denote by a^b^ and b^a^ the counterclockwise boundary arcs of dil,^ . By the Lemma IB. 21 for 
sufficiently large j the function hj is at most e{5/r) on a^b^ \ W, on the other hand having non- 
negative boundary values it is non-negative there: 

0<hj < e{6/r) on a^b^ \ W, 

and similarly 

1 - e{6/r) <hj <l on b^a^ \ W. 

Being a subsequential limit, h' also satisfies these inequalities. Sending 5 to zero, we deduce that 
h' has boundary values on ab\W and 1 on 5a \ W. 

When R goes to zero, so does r < 2R, and we see that h' has boundary values on ab and 1 on 
ba, and being bounded coincides with h. □ 

Let Q he a square with side 6L on the lattice (5Z^ and denote by 9Q a nine times bigger square. 
We will need the following continuity estimate from Verblunsky's p^7J: if a function h is discrete 
harmonic in a square 2Q, then on the square Q 

, N , const 

(23) sup |V/i| < — — sup |/i| . 

Let G(x, y) = Gqq{x, y) denote the discrete Green's function for the square 9Q, with A^G^x, y) = 
AyG(x, y) = 5x-y and vanishing on the boundary of 9Q x 9Q. By Gc we denote the discrete Green's 
function in the whole plane, normalized so that Gc{y,y) = 0. By the equation (9.6) in the paper 
|10j of McCrea and Whipple, it satisfies 



(24) Gc(x,y) = ilog^-^ + C7 + of^^V ^ - oo, 

vr 5 \\x — y\J 5 

for a specific constant C (which can be written in terms of the Euler constant). 

We will need the following integral estimate of the gradient of G in terms of G itself: 

Lemma B.4. There is a constant const independent of L such that for every y G 9Q one has 

(25) 5;|VG(x,y)|<^j;|G(x,y)|. 

xeQ x£Q 

Proof. By adjusting the constant we can assume that L is large enough. 

Suppose first that y £ 2Q. Denote by G^(-,y) the discrete harmonic function on 9Q having the 
same boundary values as Gc{-,y). We note that 

Gi;y) = Gci;y)-GU;y)- 



\x - y\ 
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By dMD on (99Q we have 

y) = Gc{; y)>l log (^^) + ^ + (z 

and so by the maximum principle the same estimate holds for G^(-,y) inside 9Q. On the other 
hand, (j24p implies that for x £ Q 

Gcix, y) < i log (2V2L) + C7 + o (^3^) • 

Combining those inequalities we infer that for x €z Q 

G(., „) < i log (2V2L) - i log (sii) + (i^) = -i log ^ + „ 

and summing over Q (recall that G is negative) we arrive at 

(26) ^\G{x,y)\>constL^. 

It follows from (fM|) that G|^(-,y) is equal on 99(5 to a constant function MogL + C up to an 
error term of ^ log (5i\/2) + o (i). Therefore by 1^ one has 

(27) J] IVG^x, y)\<Y.il (4^) + ° = 

Let £ be a lattice line through y, and i' an orthogonal line intersecting i at x and at x' and 
x" . By symmetry the whole plane Green's function G£{-,y) is monotone on the intervals [x' , x] and 
[x,x"]. So using (I24p we estimate the sum of absolute values of differences of G along this line by 

Gc{x', y) + Gc{x", y) - 2Gc(x, y) < 2 log L - 2 log + const. 



Summing this up for all lattice lines in both directions, we arrive at 

L 

(28) ^ |VGc(x,y)| < 8^ (log L - logj + const) < constL. 
xsQ j=l 

Combining (|26l) . (127^ and (I28p we prove the Lemma in the case y G 2Q: 



^ |VG| < \VGl\ + iVGcl < constL < ^ |G| . 

It remains to deal with the case y G 9(5 \ 2(5- In this case G{-,y) is discrete harmonic and 
negative in 2Q, so its values on Q are comparable by Harnack's principle to its value at the center, 
say A. Using (f23]l again, we write 



Eivgi<Et4Eigi. 



L L 

thus proving the Lemma. □ 
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Figure 6. Vertices and edges around a square. Lines corresponding to the edges 
pass through the written vectors. 

Appendix C. Unpleasant computations for Lemma 13.81 

There are several ways to prove equations (j7|8p and we present not the shortest calculation, but 
perhaps the most straight-forward one. 

We will prove ([7]), the proof of ([8]) is similar. Let u be the center of some white square. Denote 
by NW, NE, SE, SW its corner vertices, starting from the upper left and going clockwise. Recall 
that by Remark 13.31 

F{NW) - F{SE) = i {F{NE) - F{SW)) , 
so to prove ([8]) it is sufficient to show that 

AF(VF) = -\F{NE)- F{SW)\^ . 

To simplify calculations denote A := exp(— zvr/S). Denote by A^, E, S, W the centers of bordering 
edges, starting from the top and going clockwise. Assume that the line i{N) passes through a unit 
vector Q. With the chosen orientation a = 1, but we will leave a a variable to be able to compare 
results for different vertices. Then the lines £{W), i{S), i{E) pass through the vectors aX^, aA^, 
correspondingly. See Figure [H 

First evaluate increment d^wH of H across the vertex NW . Denote by A the center of the edge 
going left from the vertex NW . Recalling by definition of H we write 

dNwH = \F{W)\'-\F{A)\' 

= |Proj {F{NW),a\^)\^ - |Proj [F {NW) , a\'^)\^ 

= ^ \F{NW) + a^X^F{NW)\^ - ^ \F[NW) + a^\^F{NW)\^ 
= ^ {F{NW) + a^\'^F{NW)) {F{NW) + d?\'^F{NW)) 
- ^ {F{NW) + Q^\^F{NW)) {F{NW) + a^\^F{NW)) 
= ^ {2F{NW)F{NW) + a^X'^F^iNW) + a^X'^F^iNW)) 
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- ^ {2F{NW)F{NW) + a^\^F'^{NW) + o?\^F'^{NW)) 
= ^ ((a2 _ a6) a^\^F\NW) + - A^) a^X^F^NW)) 



^ {a^\^F\NW) + a^\^F\NW)) 



[a' 



Writing similarly increments across SW , SE, NE (with a substituted by aA^, aA^, corre- 
spondingly), we arrive at the following four equations: 



(29) 


OnwH = 


V2 
4 


{a'^X^F^{NW) + a'^X'^F^iNW)) 


(30) 


dswH = 


V2 
4 


(a'^X^F^iSW) + a'^X^F^iSW)) 


(31) 


dsEH = 


V2 
4 


{a^X^''F\SE)+a^X^''F\SE)) 


(32) 


QneH = 


V2 
4 


{a'^X^^F\NE)+a'^X^^F'^{NE)) 



Let us remark that from these equations it becomes clear that H is an appropriate discrete version 
of the primitive Im J F'^dz. Indeed, if z is one of the corners and v is the center of the square across 
that corner, denoting r/ := [v — u)/ \v — u\ we see that 

{H{v) - H{u)) = ^ {zfjF\z) - irjF\z)) 

= ^Iui{riF\z)) 

= ^Im {F\z) (v-u)), 

and therefore 

(33) 2 6 {H{v) - H{u)) = Im {F{zf {v - u)) . 

Summing the equations ()29l30l31l32p . we can write the Laplacian AH(u) in terms of the values 
of F at four neighboring vertices. But we want to reduce this further to the values of F at two 
vertices. Such reduction is possible, since by discrete analyticity projections of F{NW) on lines 
-^(A^) and i{W) coincide with those of F{NE) and F{SW) correspondingly. Using ^ we can write 
that as 

F{NW) + a'^F{NW) = F{NE) + a^F{NE) , 
F{NW) + Q^X'^FiNW) = F{SW) + a^X^FiSW) . 
Subtracting the equations multiplied by A^ and A^ correspondingly we arrive at 

(A^ - X?)F{NW) = X^F{NE) + a^X^F{NE) - >?F{SW) - a^X^F{SW) , 
where A^ — A'^ simplifies to —i\pl: 

(34) F{NW) = ^ {X'^F{NE) + a^X^F{NE) - X^F{SW) - a^X^F{SW)) . 

v2 
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Writing similarly for F{SE) (with aA^ substituted for a and NE and SW interchanged) we con- 
clude that 

(35) F{SE) = {X^F{SW) + a'^X^^F{SW) - X'^F{NE) - a^X^^F{NE)) , 

v2 

Now we can sum equations (l29l31l3UI32p . substituting (|34I35|) for values of F{NW) and F{SE): 
AH{w) = dNwH + dswH + dssH + OneH 

= ^ (^a^X^ {X^F{NE) + a^X'^F{NE) - X'^F{SW) - a^X'^F{SW)) ^ 

+ 6?X^ (^-^ {X^F{NE) + a^X^F{NE) - X^F{SW) - a^X^F{SW))^ 
+ a^A^o _!= {X^F{SW) + a2AiOF(5i^) _ X'^F{NE) - a'^X^OF{NE)) ^ 
+ a^A^o (^-^ {X^F{SW) + a'^X^°F{SW) - X^F{NE) - a^X^^F{NE))^ 
+ a^X^F^{SW) + a'^X^F^iSW) + a^X^^P{NE) + a^X^^F'^{NE)^ 

When we plug in a = 1 and recall that A = exp(— z7r/8), in particular A*^ = —1, there are many 
cancelations in the right hand side, which eventually simplifies: 

■ ■ ■ " X (^^ ( "7! + X^F{NE) - X^F{SW) - X^F{SW)) ^ 

+ P (^-^ {X^F{NE) + X'^F{NE) - X'^F{SW) - X^F{SW))^ 

- A2 --^ {X'^F{SW) - X'^F{SW) - X'^F{NE) + X^F{NE)) ^ 

- A2 {X^F{SW) - X'^F{SW) - X'^F{NE) + X'^F{NE))^ 
+ X^F^iSW) + X^F^iSW) - X^F^iNE) - X^F^{NE) 



V2f , .3p.™ T^.™\2 



X V 2 + ^^i^E) - ^ F{SW) - XF{SW) 

- i (^XF{NE) + XF{NE) - X^F{SW) - XFiSW)^ 
+ ^ (lF{SW) - XF{SW) - X:'^F{NE) + XF{NE)'^ ^ 
+ l(XF{SW) - XF{SW) - X-^FiNE) + XF{NE) 

^1 V 
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= ^ (^F\NE) (-A^ -X^ + X^ + X^- 2X') + F\NE) - + + - 2A6) 

+ F'^{SW) (-A^ - A'^ + A^ + A^ + 2A6) + F'^{SW) (-A^ - A^ + A^ + A^ + 2A^) 
+ 2F{NE)F{SW) (A^ + A^ - A^ - A^) + 2F{NE)F{SW) (A^ + A^ - A^ - A^) 
+ 2F{NE)F{NE) (-A^ - A^ - A^ - A^) + 2F{SW)F{SW) (-A^ - P - A^ - A^) 

+ 2F{NE)F{SW) (A^ + A^ + A^ + A^) + 2F{SW)F{NE) (A^ + A^ + A^ + A^) j 

/2 

= -^4V2{F{NE) - F{SW)) {F{NE) - F{SW)) = - \F{NE) - F{SW)f . 
8 

This finishes the proof of Lemma. 
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